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Abstract 

FeynHiggs is a Fortran code for the calculation of the masses of the neutral CP- 
even Higgs bosons in the MSSM up to two-loop order. It is based on the complete 
diagrammatic on-shell results at the one-loop level, the leading diagrammatic two-loop 
QCD contributions and further improvements taking into account leading electroweak 
two-loop and leading higher-order QCD corrections. The Higgs-boson masses are cal- 
culated as functions of the MSSM parameters for general mixing in the scalar top sector 
and arbitrary choices of the parameters in the Higgs sector of the model. 
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1 Introduction 



A direct and very stringent test of Supersymmetry (SUSY) is provided by the search for 
the lightest Higgs boson, since the prediction of a relatively light Higgs boson is common to 
all supersymmetric models whose couplings remain in the perturbative regime up to a very 
high energy scale PJ. A precise prediction for the mass of the lightest Higgs boson, m^, in 
terms of the relevant SUSY parameters is crucial in order to determine the discovery and 
exclusion potential of LEP2 and the upgraded Tevatron. If the Higgs boson exists, it will be 
accessible at the LHC and future linear colliders, where then a high-precision measurement 
of the mass of this particle will become feasible. A precise knowledge of the mass of the 
heavier CP-even Higgs boson, m#, is important for resolving the mass splitting between the 
CP-even and -odd Higgs-boson masses. 

In the Minimal Supersymmetric Standard Model (MSSM) [Q] the mass of the lightest 
Higgs boson is restricted at the tree level to be smaller than the Z-boson mass. This bound, 
however, is strongly affected by the inclusion of radiative corrections, which yield an upper 
bound of about 130 GeV J|, [|, ^ || 0, ||, [|. Results beyond one-loop order have been 
obtained in several approaches: a Feynman diagrammatic calculation of the leading QCD 
corrections has been performed || [7]]. Also renormalization group (RG) methods have been 
applied in order to obtain leading logarithmic higher-order contributions [§] . Furthermore the 
leading two-loop QCD corrections have been calculated in the effective potential method ||. 
All these calculations show that the corrections beyond one-loop order lead to a sizable 
decrease of m h of up to 20 GeV. 

Concerning the calculation of the lighter and heavier neutral CP-even Higgs bosons two 
different kinds of computer codes have been used for phenomenological analyses so far: they 
are either based on the RG improved one-loop effective potential approach || or on the 
one-loop diagrammatic on-shell calculation [§|, |J. These approaches differ by 0(10 GeV). 

Here we present a new Fortran program named FeynHiggs, which is based on the results of 
the Feynman-diagrammatic calculations up to (D(aa s ) given in Refs. 0, ||, 0]. It includes the 
complete diagrammatic on-shell results at the one-loop level, the leading diagrammatic two- 
loop QCD contributions and further improvements taking into account leading electroweak 
two-loop and leading higher-order QCD corrections. The calculation of mh and win is 
performed for arbitrary values of parameters in the t — t-sector and the Higgs sector of 
the MSSM. As a subroutine, FeynHiggs can be linked to other programs thus incorporating 
in an easy way a precise prediction for rrih and m#. In addition the program provides as 
an option the calculation of the SUSY contribution to the p-parameter in 0(aa s ), based 
on Refs. [JIT], [□]]. In this way experimentally disfavored combinations of squark masses can 
automatically be excluded. 

The paper is organized as follows: in section 2 we specify our notations and give a brief 
outline of the calculation of nth and m#. In section 3 the program FeynHiggs is described in 
detail. Examples of how to use FeynHiggs are shown in section 4. The conclusions are given 
in section 5. 
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2 The calculational basis 



2.1 The top-squark sector of the MSSM 

In order to fix the notation we shortly list our conventions for the MSSM scalar top sector: 
the mass matrix in the basis of the current eigenstates ti and is given by 

u2 _( Mf L + ml + cos2(3 (I- hw)M 2 z m t M t LR \ m 

Mi ~ { m t M t LR M? r +m 2 t + \ cos 2/3 s 2 w M 2 z J ' [L) 

where 

m t M t LR = m t (A t - // cot (3 ) . (2) 

Diagonalizing the t-mass matrix yields the mass eigenvalues , m^ 2 and the t mixing an- 
gle 8t, which relates the current eigenstates to the mass eigenstates: 



/ ti \ _ ( cos^- sin#£ \ / ti 
\ i 2 J V - sin 9- t cos d~ t J \t R 

2.2 Calculation of the Higgs-boson masses 



(3) 



Contrary to the Standard Model (SM), in the MSSM two Higgs doublets are required. The 
Higgs potential JT2 



V = mlH 1 H 1 + mlH 2 H 2 -m 2 12 {e ab H'tH* + h.c.) 

+ {HA - H 2 H 2 f + £ \H X H 2 \ 2 (4) 

contains mi, m 2 , m\ 2 as soft SUSY breaking parameters; g, g' are the SU (2) and U{1) gauge 
couplings, and e\ 2 = —1. 

The doublet fields Hi and H 2 are decomposed in the following way: 



rr (Hi) _ ( Vl + {(j) + ix )/V 2\ 



H2 ~ (ff!)"(«2 + (^ + <^)/v^)- (5) 

The potential (f|) can be described with the help of two independent parameters (besides 
g and g')\ tan/3 = v 2 /v\ and M\ = — m^ 2 (tan/5 + cot (3 ), where Ma is the mass of the 
CP-odd A boson. 

In order to obtain the CP-even neutral mass eigenstates, the rotation 



H° \ ( cos a sin a \ / <fi\ 
h° J I — sin a cos a I \ 4>2 



(6) 



is performed, where the mixing angle a is given in terms of tan (3 and Ma as follows: 



Ml + Ml 
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tan 2a = tan 2(3^ -§, < a < 0. (7) 



MX~M 2 z 2 



At tree level the mass matrix of the neutral CP-even Higgs bosons in the 0i — 02 basis 
can be expressed in terms of Mz and Ma as follows: 



M 2 ' tree - ( m <t>i m 0i<fe 1 
m m ggs — I 2 m 2 

( M\ sin 2 /? + Mf cos 2 -(M| + M|) sin /3 cos /3 \ 

^ -(Ml + M|) sin/? cos/? M| cos 2 /? + M| sin 2 (3 )' W 

which by diagonalization according to eq. (^) yields the tree-level Higgs-boson masses 

M 2,tme ( m H,tree | (q\ 

Mffisgs m 2 • ® 



In the Feynman-diagrammatic approach the Higgs-boson masses in higher orders are 
derived by finding the poles of the h — //-propagator matrix whose inverse reads: 



, x x-i ■( Q 2 - m H. tree + ^H(q 2 ) £hH(q 2 ) \ M , 



The poles are then obtained by solving the equation 

- Ktree + EfcfaW " <,tree + M^)) - (W<? 2 )) 2 = . (11) 

In the following the fiw denote the one-loop (i = 1) and the two-loop (i = 2) contributions 
to the renormalized self-energies. 

In FeynHiggs the one-loop results for the Higgs-boson self-energies E^(g 2 ) are calculated 
according to Ref. ||. They contain the full one-loop contribution obtained via an explicit 
Feynman-diagrammatic calculation in the on-shell scheme. Here the gaugino parameters Mi 
and M(= M 2 ) enter in the neutralino mass matrix. Mi is fixed via the GUT relation 

Mi = l&M, (12) 
"J &w 

whereas M is kept free input parameter. 

The two-loop results for the Higgs-boson self-energies are taken from Refs. || 0] . The 
leading two-loop corrections have been obtained by calculating the 0(aa s ) contribution of 
the t— t-sector to the renormalized Higgs-boson self-energies at zero external momentum from 
the Yukawa part of the theory. These two-loop QCD corrections are expected to constitute 
the most sizable part of the full set of two-loop corrections. In Refs. [pi U\ the self-energies 
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have been computed first in the <\>\ — 02 basis and afterwards rotated into the h — H basis 
according to eq. (f|): 



M2) 
-'H 

-'h 

s(2) 

J hH 



cos a + sin a + 2 sin a cos a 
sin a h + cos a L^ 2 — 2 sm a cos a 



— sm a cos a 



(Eg - Eg) + (cos 2 a -sin 2 «)S 



(2) 



(13) 



Thus for our results up to the two-loop level the matrix (|T0D contains the renormalized 
Higgs-boson self-energies 



± s (q 2 )=±^(q 2 )+±f\0), s = h,H,hH, 



(14) 



where the momentum dependence is neglected only in the two-loop contribution. The cal- 
culation is performed for arbitrary parameters of the Higgs and the scalar top sector and for 
arbitrary gluino mass rrig. Thus the accuracy of the calculation does not depend on how the 
parameters of the t sector mj , mx and 8z are chosen. 



In order to take into account the leading electroweak two-loop contribution to the mass of 
the lightest Higgs boson, we have implemented the leading Yukawa correction of 0(G 2 F m^), 
which gives a sizable contribution only for m^. The formula is taken over from the result 



obtained by renormalization group methods. It reads [13 



Am 2 



with 




sin 2 26; 



1 2 

mi —mi 
f \ „ fl sin 2 2^-log 



log 



2m\ 

'm~m~ ' 
mf 



m~ + m~ 

Hf H°S 

mi — mi 





(15) 



(16) 
(17) 



The second step of refinement incorporated into FeynHiggs concerns leading QCD cor- 
rections beyond two-loop order. They are taken into account by using the MS top mass 



m t = m t [m t ) 



m t 



nF a »( m *) 



for the two-loop contributions instead of the pole mass, mt = 175 GeV. 

The results implemented in FeynHiggs have been compared to the calculations using RG 
methods. Good agreement has been found in the case of no mixing in the t sector, i.e. 
M t = GeV, whereas sizable deviations can occur when mixing in the t-sector is taken 
into account. This has been discussed in detail in Ref. H. 



2.3 Calculation of the /^-parameter 



We have also implemented the calculation of the MSSM contributions to Ap |fT0| , [LI]. Here 
the corrections arising from i/6-loops up to 0(aa s ) have been taken into account. The result 
is valid for arbitrary parameters in the t- and 6-sector, also taking into account the mixing 
in the 6-sector which can have a non- negligible effect in the large tan/3 scenario |TT|| . 

The two-loop result is separated into the pure gluon-exchange contribution, which can 
be expressed by a very compact formula that allows a very fast evaluation, and the pure 
gluino-exchange contribution, which is given by a rather lengthy expression. The latter 
correction goes to zero with increasing gluino mass and can thus be discarded for a heavy 
gluinoQ The p-parameter can be used as an additional constraint (besides the experimental 
bounds) on the squark masses. A value of Ap outside the experimentally preferred region of 
Ap SUSY « 1CT 3 [P| indicates experimentally disfavored t- and 6-masses. 



3 The Fortran program FeynHiggs 
3.1 The main structure 

The complete program FeynHiggs consists of about 50.000 lines Fortran code, where the 
main part belongs to the formulas for the renormalized Higgs self-energies at the two-loop 
level and the gluino contribution to Ap. The executable file fills about 4 MB disk space. 
The calculation for one set of parameters, including the Ap constraint, with the highest 
accuracy takes less than 2 seconds on a third-generation Alpha 21164 microprocessor (300 
MHz processing speed). 

There exists a Home page for FeynHiggs: 
bvttp : / /www-itp . physik . uni-karlsruhe . de/f eynhiggs . 

Here a uu-encoded version as well as an ASCII version is available, together with a short 
instruction, information about bug fixes etc. 

FeynHiggs consists of several subprograms which are listed in Table |l|. We now describe 
the different subprograms in detail. 

• FeynHiggs. f is the front-end for the whole program. Here all variables are set, all 
options are chosen, the subprogram FeynHiggsSub.f is called, and the results for 
the Higgs masses are printed out. For an easy use of FeynHiggs this front-end can be 
manipulated at will, the rest of the program need not to be changed. 

• FeynHiggsSub.f: Here the actual calculation is carried out. The various results for 
the renormalized Higgs-boson self-energies ( |14"D are put together. Eq. (|TT|) is solved 
numerically, the refinement terms (the leading two-loop Yukawa contribution and the 
leading QCD corrections beyond two-loop order) are incorporated. 

1 An additional contribution to Ap, arising from a shift in the squark masses when the soft SUSY breaking 
parameters are used as input (due to the SU(2) invariance of these parameters in the squark sector), is not 
implemented in FeynHiggs. This correction is described in detail in Ref . [fill . 
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subprogram 


function of the program 


FeynHiggs.f 


front-end 


FeynHiggsSub.f 


main part of FeynHiggs : calculation of m^mu 


Hhmasssr2.f 


one-loop self-energies 


varcom.h 


definition of variables 


bc.f 


one-loop functions 


lamspen.f 


further mathematical functions 


def2.f 


definitions for the MSSM parameters 


Plsecode.f 


Fortran code for S^(0) 


Plsesum.f 


putting together the different parts of £^(0) 


Plsevar.f 


definition of variables for S^(0) 


P2se[code,sum,var].f 


same for E^(0) 


PlP2se[code,sum,var].f 


same for E^ 2 (0) 


Plsetl[code,sum,var] .f 


same for EjJ(O) 


P 2 set 1 [co de , sum, var] . f 


same for sJJ(O) 


PlP2setl[code,sum,var].f 


same for £^ 2 (0) 


delrhosub.f 


main part for the calculation of Ap 


delrhoGluino[code,sum,var] .f 


subroutines for the gluino-exchange contribution to Ap 



Table 1: The subprograms of FeynHiggs. 



• Hhmasssr2.f: The results for the complete one-loop self-energies are evaluated. 

• varcom.h: The variables needed for the files *code.f are defined and grouped in 
common blocks. 

• bc.f, lamspen.f contain mathematical functions needed for the one- and two-loop 
self-energies. 

• def2.f: The masses of the stops and sbottoms are calculated from the parameters in 
the squark mass matrices. The opposite way is also possible. 

• Plsetlcode.f contains the complete code for the two-loop contribution to the renor- 
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malized <f>\ self-energy in 0(aa s ). This code was first calculated with the help of 
Mathematica packages (see || |7j) and was afterwards transformed automatically into 
this Fortran code. The complete code has been split up into 13 smaller functions in 
order not to exceed the maximal number of continuation lines of the Fortran compiler. 

• Plsetlsum.f: The sum of 13 functions for the renormalized <pi self-energy at the 
two- loop level, contained in Plsetlcode.f, is put together. 

• Plsetlvar.f contains the variable definition for the 13 subfunctions described above. 



• The above three files exist in an analogous way for the other renormalized Higgs-boson 
self-energies at the one- and at the two-loop level. The one-loop self-energies in these 
files are given in the approximation explained in detail in [0, |7J and are only needed 
internally for consistency checks. For the real calculation of my, and run the complete 
one-loop self-energies are calculated in Hhmasssr2.f. 

• delrhosub.f is needed for the calculation of the MSSM contribution to Ap. In this 
program the evaluation of the leading one-loop contribution is performed. In addition 
also the gluon-exchange correction in 0{aa s ) is calculated. The 0(aa s ) contribu- 
tion to Ap is completed by including the gluino-exchange contribution by calling the 
subprogram delrhoGluinosum.f. 

• delrhoGluino [code, sum, var].f: The code for the gluino exchange contribution to 
Ap is implemented in the same way as it is described for the renormalized Higgs-boson 
self-energies. 

3.2 Options and setting of variables 

FeynHiggs can be run in several ways, determined by the choice of several options. 



Depth of calculation allows to choose to what extent the refinements described in 



section 2.2 should be applied. 



Selection of input parameters: One can either use the physical parameters of the 
t mass matrix (m^ , m^ 2 and or the soft SUSY breaking parameters (M$ L , M^ R and 
M^ R ) as input parameters. 

m t in the t mass matrix: If the soft SUSY breaking parameters are used as input 
parameters, one can choose whether the top pole mass, m t , or the running top mass, 
Tn t , should be used in the t mass matrix in order to determine the masses of the 
eigenstates m^rn^ and 9^. 

Limit for Ap SUSY : One can specify the maximally allowed value for the MSSM con- 
tribution to Ap. If Ap SUSY exceeds this limit a warning is printed out. 
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• Selection for the one-loop accuracy: Before the calculation of the Higgs masses 
starts, one has to specify to what accuracy the one-loop renormalized self-energies 
should be evaluated. One can either take into account the top sector only, one can 
choose to use the top and the bottom sector, or one can select the option that the 
complete MSSM should be taken into account. 



input for FeynHiggs 


• • it TV If f~T TV If 

expression in the MSSM 


internal expr. in FeynHiggs 


tan (beta) 


tan (3 


ttb 


Msusy_top_L 




msusytl 


Msusy_top_R 


M ~tn 


msusytr 




M h 


msusybl 






msusybr 


MtLR 


M t LR 


mtlr 


MSt2 


m h 


mst2 


delmst 


Am~ t = m h - m h 


delmst 


sin (theta_stop) 


sin 9i 


stt 




sin 61 


stb 


Mtop or mt 


m t 


mmt 




m b 


mbb 


Mgluino 


m- g 


mgl 


Mue 


n 


mmue 


M 


M 


mmm 


MA 


M A 


mma 



Table 2: The meaning of the different MSSM variables which have to be entered into Feyn- 
Higgs. 



In FeynHiggs. f the Standard Model (SM) variables are set to their present experimental 
values. New values for these variables can be implemented by the user into the code of the 
file FeynHiggs. f easily. The MSSM variables can be chosen by the user at will. In addition, 
since the dependence of the Higgs masses on the top mass is very strong, also mt can be 
chosen at will. In Table |2] the meaning of the different parameters FeynHiggs asks for is 
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explained^. All these variables are transferred to the different subprograms by common 
blocks. 



4 How to use FeynHiggs 

In this section we will give two examples of how FeynHiggs is used. As stated before, the 
front-end FeynHiggs. f can be manipulated by the user at will, whereas the subprogram 
FeynHiggsSub.f and all the other subprograms should not be changed. Concerning a 
modificaton of the front-end one has to note that all variables have to be defined in the 
front-end. The Higgs masses are then obtained by calling the subroutine via 

call f eynhiggssub(mhl, mh2, mhl2, mh22) , 

where mhl and mh2 are the one-loop corrected values for rrih and m#, respectively. mhl2 and 
mh22 contain the values for the two-loop corrected Higgs masses, including the refinement 
terms as it has been specified in the options. 

In the following two examples are given, how the application of FeynHiggs looks like, 
using the given front-end of the currently distributed version. The user's input is given in 
bold face letters. 



2 Some MSSM variables exist internally in FeynHiggs, but are not input parameters. These variables 
have no entry in the left column of Table 
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4.1 Example 1 



>FeynHiggs.exe 
. . . Introduction . . . 



depth of calculation ? 

1: full 1-loop + 2-loop QCD 

2: same as 1, but in addition with 

mt = 166.5 at 2-loop 
3: same as 2, but in addition with 

Yukawa term added for light Higgs 

1 

Select input: 

1 : Msusy, MtLR, . . . 

2: MSt2, delmst, stt, ... 

2 

Limit for Delta rho = 1.3 * 10" -3 ? (0 = ok) 




tan (beta) = ? 
1.6 

MSt2 = ? 
400 

delmst = ? 
100 

sin(theta_stop) = ? 

(0: stt = // 1: stt = sin(-Pi/4)) 

1 

Mtop = 175 ? (0 = ok) 
173.8 

Mgluino = 500 ? (0 = ok) 
300 

Mue = -200 ? (0 = ok) 
100 

M = ? (0: M = 400 // 1: M = Msusy) 
1 
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MA = ? 
500 

Selection: 1 = top only, 2 = top/bottom only, 3 = all 
3 

Your parameters : 

tb, Msusy(top-lef t) , Msusy(top-right) , MtLR 
MT, Mgl, Mue, M, MA 

1.600000 309.2813 308.0858 200.0000 

173.8000 300.0000 100.0000 309.2813 500.0000 



The results: light Higgs 



heavy Higgs 



mh-tree 

mh-lloop 

mh-21oop 



39.42879 
69.53253 
64.10773 



506.7153 
508.7314 
508.3541 



Delta rho 1-loop 
Delta rho 2-loop (gluon) 
Delta rho 2-loop (gluino) 
Delta rho total 



5 . 909909368803866E-004 
6 . 177715322384914E-005 
4 . 755834893556049E-006 
6 . 575239249977918E-004 



tan (beta) = ? 

In this example the physical parameters in the i-sector have been chosen as input pa- 
rameters, no refinement term has been included. The selection for M sets M = Mi , where 
Mi L is calculated from , m 4 ~ 2 and sin 9$, see eq. ([[J) . 
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4.2 Example 2 



>FeynHiggs.exe 
. . . Introduction . . . 



depth of calculation ? 

1: full 1-loop + 2-loop QCD 

2: same as 1, but in addition with 

mt = 166.5 at 2-loop 
3: same as 2, but in addition with 

Yukawa term added for light Higgs 

3 

mt in the stop mass matrix at 2-loop ? 
1 : mt = pole mass 
2: mt = running mass 
2 

Select input : 

1 : Msusy, MtLR, . . . 

2: MSt2, delmst, stt, ... 

1 

Limit for Delta rho = 1.3 * 10" -3 ? (0 = ok) 
0.00001 



tan (beta) = ? 
20 

Msusy_top_L = ? 
1000 

Msusy_top_R = ? (0: Msusy_top_R = Msusy_top_L) 


Mtop = 175 ? (0 = ok) 


Mgluino = 500 ? (0 = ok) 


Mue = -200 ? (0 = ok) 
-100 
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M = ? (0: M = 400 // 1: M = Msusy) 

300 

MA = ? 

400 

MtLR = ? 
1000 

Selection: 1 = top only, 2 = top/bottom only, 3 = all 
3 

Your parameters : 

tb, Msusy (top-left) , Msusy(top-right) , MtLR 
MT, Mgl, Mue, M, MA 

20.00000 1000.000 1000.000 1000.000 

175.0000 500.0000 -100.0000 300.0000 400.0000 



The results: 

mh-tree : 
mh-lloop : 
mh-21oop : 



light Higgs 

90.70748 
133.6753 
118.8657 



heavy Higgs 

400.1090 
400.1799 
400.1770 



using running mt for two-loop contribution: 167.338636349986 
. . . also for mt in Stop mass matrix 
mh-21oop: 120.8774 400.1780 



using running mt for two-loop contribution: 167.338636349986 

. . . also for mt in Stop mass matrix 

adding Yukawa term for light Higgs 

. . . also with running mt in stop mass matrix 

mh-21oop: 122.2504 400.1780 



WARNING: Delta rho > experimental limit 
Delta rho 1-loop : 3 . 224156596235517E-005 

Delta rho 2-loop (gluon) : 3.475299800114144E-006 
Delta rho 2-loop (gluino) : 2 . 896993903992095E-005 
Delta rho total : 6.468680480239026E-005 



tan (beta) = ? 
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If one selects to include all refinement terms, the front-end automatically calculates 
the Higgs masses in all three steps of accuracy. The running top mass (in this example 
167 . 338636349986 GeV) is used also in the t mass matrix. In this example only a very 
small SUSY contribution to Ap ia allowed and the value of Ap Y for the chosen parameters 
exceeds the above specified maximal value (a warning is printed out.) 

5 Conclusions 

FeynHiggs is a Fortran code for the calculation of the masses of the neutral CP-even Higgs 
bosons of the MSSM. It is based on results which have been obtained using the Feynman- 
diagrammatic approach. The results consist of the full one-loop and the leading two-loop 
QCD corrections. Two further steps of refinements have been implemented. The Fortran 
code for the one-loop contribution has been taken over from a program written by A. Dabel- 
stein. The Fortran code for the two-loop correction has been generated automatically from 
a Mathematica result. 

The program is available via the WWW page 
http : / /www-itp . physik . uni-karlsruhe . de/f eynhiggs . 

The code consists of a front-end and a subroutine. The front-end can be manipulated 
at the user's will. The different parts of the code have been described in detail, and the 
meaning of the variables used in the code has been explained. We have given two examples 
of how to use FeynHiggs. 

The subroutine is self-contained and can be linked as an independent part to other 
existing programs^ thus providing a precise prediction for rrih and m#, which can then be 
used for further computations. 
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